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PREFACE 


As part of the Ministry of the Environment's contribution 
to the International Hydrological Decade, the River Basin 
Research Section is carrying out studies of water resources 
and physical conditions in five basins in Southern Ontario. 
The basins were selected as being representative of common 
type areas in the province and the hydrologic studies being 
undertaken are designed to provide a better understanding of 
most aspects of the water balance in these areas. 


This paper outlines a proposed approach towards the 
development of mathematical models that will describe the 
complex interactions of the various components of the hydrologic 
cycle in each basin. Such models will lead to a better under- 
standing of the actual hydrologic processes that exist and will 
be useful in the synthesis of hydrologic characteristics in 
ungauged areas thereby aiding in the solution of water resources 
problems. 


It is anticipated that the compilation of pertinent 
information and approach outlined in this paper may be helpful 


to others engaged in the development of mathematical hydrologic 
models. 


K. EB. Symons, Director, 
Water Quantity Management Branch. 


Toronto, November 1, 1972. 
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AN APPROACH TO MATHEMATICAL MODELLING OF 
MINISTRY OF THE ENVIRONMENT 


IHD REPRESENTATIVE BASINS 


INTRODUCTION 
According to the 'Guidelines for Research Basin Studies', 

as compiled in 1966 by the Canadian National Committee for the 
International Hydrological Decade, 'representative basins are 
areas set aside to permit study of the existing conditions of 
a region'. One of the objectives of research in these areas 
is: 

'To provide data, on all parameters in 

the hydrologic cycle, so that statistical 

or physical models may be developed for 

major climatic and geologic areas. * 
Principles derived from these models will 


be used to synthesize the hydrological 


characteristics of ungauged areas.' 


An aim of hydrologic research is to interpret collected 
data for use in the solution of water management problems. 
All drainage basin studies should therefore be planned to 
further the understanding of the hydrologic cycle by using 


knowledge of the physical and mathematical relationships 


* See Appendix I for location maps of basins representative 
of major climatic and geologic areas in southern Ontario, 
under study by the River Basin Research Section. 
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ok 
between the various components. The development of ‘models' 
from these relationships should then provide the ability to 
accurately reconstruct past and present conditions and to predict 
future events, relating to hydrological phenomena. These models 
can therefore be useful aids to solving water resources problems 
dealing with evaluation, design, management and regulation. 
Hydrologic Models 

Runoff from homogeneous drainage areas and Pores variant flows 
in open channels are affected by factors which are generally 
known and can be defined explicitly. These factors are difficult 
to evaluate in a basin-wide study where data are usually scarce 
and vary both in time and space. However, the analysis of a 
rainfall-runoff process in a complex drainage basin can be made 
feasible by the use of models based on concepts of the system. 
Basin models may be classified into three general categories: 

a) physical models, 

b) analog models, 

c) mathematical models. 
In the selection of a type of model, consideration must be given 
to cost of development, operational problems, accuracy, inter- 
pretation of results and subsequent model flexibility. 

Physical models are used to represent an idealized physical 
process or environment and are designed for a specific purpose, 
at a specific location. Physical models are employed in 
hydrology to simulate such processes as infiltration, ground- 


water flow and rainfall-runoff phenomena. With minor modifications, 


aay. oes ate path itt 
; CATS) eT 
oo, 
; 50 ug ice, ol pie vik i ia sin iis 
, : _ a yo : : . 
- « deltee a4 oF byis dap} keno “nse: *q feeg wae 
hae Ds 7 ear ‘ 
vt eere rec Lan re :. a 
Y U4 shaman bier .sasmonbidg Leake oat oF cuales 2 
¥ re aye | dosyoacizuqes: taise yprty.os’ os abts sacilnieed a i 
“sh xoidefopes 6 4ns! ine pce olvaiitsve tk 


io tube A okt 


aol? InaTier-srts ns aire venmant move: 


' s 
/ , = te 
La = 
; - i t ya? } 
ian F; 
o - b p | a 
: 
I nica 1 
s) OF ST t7 
’ yy 
af 
| 
a5 
+ 
\ 
" 


, Bae \ pit vite Woes MNES He 
=) [an eT LAL f es MOoltsool] ora 


Uns & 


,enol2 Se zon ar 


these models can often be used for extensive periods of time, 
justifying their cost of construction and operation through the 
importance of such economic essentials as water-supply, flood 
control and navigation. 

Analog models are commonly used to a varying degree in 
solving ground-water flow problems. The principle involved in 
this type of modelling utilizes the analogy between the flow of 
electricity in an electronic circuit and the flow of water in 
an aquifer. Analog models, like physical models, are usually 
designed for a specific location, with the costs being justified 
primarily on a single#purpose basis. 

Mathematical models allow natural hydrologic phenomena to 
be investigated as systems. The systems approach enables the 
development of quantitative relationships between the various 
hydrologic variables which are usually expressed through 
mathematical formulations (models) for simulation purposes. If 
the hydrologic variables are considered to follow physical laws 
rather than laws of probability, the model is described as 
parametric or deterministic. If laws of probability are considered 
in characterising the hydrologic variables, the model is described 


as stochastic or probabilistic. 


Parametric Modelling 


In preliminary hydrologic studies, parametric models of 
surface runoff usually related runoff to rainfall and drainage 
area, and took the form of a purely empirical relationship or 


of a design procedure following the rational method. With the 
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increase of knowledge in the mechanics of physical hydrologic 
processes, modelling techniques become more complex; hence, 
with a more comprehensive approach, parametric modelling can 
be treated on the basis of component modelling and integrated 
system modelling. 

(a) Component Modelling: 

The overall hydrologic cycle can be divided into several 
component processes such as precipitation, snowmelt, soil 
moisture flow, evapotranspiration, direct surface runoff, 
baseflow and total runoff. Each component or combination of 
components, may be treated as a separate model. For each model, 
empirical approximations have been developed to simulate the 
processes controlling the component(s). However, by gaining 
an improved understanding of the physical laws governing a 
component, improvements on the initial empirical models have 
been achieved whereby the governing concepts underlying the 
physical situations can be more closely approximated. 

(b) Integrated System Modelling: 

Through the development of better component models, the 
Overall modelling of the integrated or total system of a watershed 
has become more complex. Through the use of electronic and 
digital techniques, the limitations on the complexity are 
alleviated and the individual component models can be linked 


to form a comprehensive, integrated system model. 
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Stochastic Modelling 


The approach in stochastic modelling is to present a 
statistical simulation of a measured response of a system. 
In order to use effectively a stochastic model, a probability 
distribution for the input must be determined or assumed. In 
past hydrologic studies, considerable effort was made to 
determine the probabilistic properties of hydrologic data, 
particularly those of streamflow and rainfall. 

One of the characteristic features of the stochastic 
method is that the statistical measures of the hydrologic 
variables used in the calibration of the model enable probability 


limits of accuracy to be placed on the simulated values. 
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OBJECTIVES 

The general purpose of this paper is to present the 
conceptualization in hydrologic modelling relative to present 
and projected applications to hydrologic investigations in 
Ministry of the Environment - IHD representative basins. The 
schematizations developed will serve as flexible frameworks 
for a subsequent detailed approach and a step-by-step procedure 
in the development of applicable hydrologic models. 

The encompassing objective of the hydrologic simulation 
studies is to formulate models which possess detailed and 
specific features which are characteristic of a representative 
Study area. The models that are formulated should possess 
equally sufficient generality for extrapolation to other study 
areas or drainage basins within the larger study region. 
Detailed information, hydrologic and geologic data, compiled 
for the selected representative basins, are to be used in the 
optimization of the model parameters and the overall model 


calibration. 
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HYDROLOGIC MODELLING 

The review of methods, approaches and procedures in hydro- 
logic modelling, which are to be used as bases and guidelines 
in the development of working models for each of the Ministry 
of the Environment - IHD representative drainage basins, is 
summarized in the appendices. These summaries are outlined in 
a series of figures as model schematizations. The following 
text serves to introduce briefly the content of each Appendix. 
Accompanying each Appendix is a further elaboration on the text 
which serves to explain in more detail the schemes of the 
respective figures. 

These hydrologic modelling procedures are outlined on the 
basis of: 

1. Basic Conceptualization 

exemplified by, 

rk) A schematization for functional and inter-related 
hydrologic processes. 

ii) A schematization for the superposition of basin 
instrumentation on the functional and inter- 
related hydrologic processes. 

2. Component Models 


exemplified by, 


1) A schematization for hydrologic component models 


for data generation and simulation. 
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3. An Integrated Hydrologic System Model 
exemplified by, 
i) A schematization for an integrated (comprehensive) 
hydrologic system model. 
4. A Specific Model Approach 
exemplified by, 
i) A definition scheme for a parametric, linear 
regression model. 
ii) Expanded definition scheme for a parametric, 
linear regression model. 
iii) Flow chart for a parametric, linear regression 
model. 
5. A Brief on Stochastic Models 
exemplified by, 
i) A schematization for a type of stochastic 


hydrologic system model. 


Basic Conceptualization 


The basic conventional concepts of the hydrologic cycle 
were used to formulate the schematization of Figure 1 of 
Appendix II. The functional and inter-related hydrologic 
processes are outlined by means of a flow scheme showing 
basically the paths precipitation, as input, would follow 
through the schematization to the various outputs. 

The standard, measured meteorological parameters, precip- 


itation, P, temperature, T, radiation, R, -- (longwave radiation, 
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Ri, shortwave radiation, Ror Or net radiation, Rs duration 
of sunshine, S, humidity, H, wind, W, and Class-A pan eva- 
poration, are shown as linkages between the source and 
the processes for which they are a cause or effect, or on 
which they have an influence directly or indirectly. 

The precipitation process resulting in rain and/or snow 
en route to supply the surface storages, Liquid Amount and 
Snowpack, is affected by the physical process of Interception 
which results in a reduction in the total amount available to 
these surface storages. In the process, an amount is trapped 
by the vegetal canopy as Interception Storage. The surface 
storage, Snowpack (snow), is subsequently affected by certain 
physical interacting processes (heat exchange and melting), 
resulting in equivalent liquid amounts. The momentary storage 
at the surface as rainfall amounts and the melt drained from 
the snowpack are the initializing sources of input to the 
Soil Profile Storage, a subsurface storage, which is affected 
by the physical process of Infiltration. 

That portion of the initializing input which is in excess 
of the infiltration requirement, goes to initialize the surface 
physical process, Surface Runoff (overland flow) by way of 
satisfying the Depression and Detention Storages. Eventually, 
amounts from the Depression and Detention Storages, infiltrate 


to the Soil Profile Storage. 


Re sane & sexs wie dak 


rs 
— 4 YLtneEt speertiat oes ee 
By 


f, 
<7 
iL 


Pekin sites itgeeda mative’ ay yloqva £ 


. gin ieee So egenorg fetieyvds of¢ vd hessottes a 
F 
os ai waitin ‘nye, fated alt i aoitoubez ss ab af, 
7 barguys: al 2nyerce 1% teens af? al 2eosrote bad 
| Seras oA se TIR Noltgesrstol us vaonss Lidegale: 
A = 
Martens Wo Satsedte yl tirecpesedve (wore) Jb sawone = E 
,onts Ler 286 Genttinxs tour) nsanot fo -pits *"s2etnt take 
S6erutes Yietnwnat ost e-em ofepil inolaviuss 
Got? Dberiat6 gf as 2 Dae Sfanwwm iistotes Bs 
rt) «ag Pos. 70 Annas pal fPisittint: By ¢ 
BOsNHE2e SL HatilW .opstuse ‘sas 37 uz | £ ,spALOSE 
i f « I ‘, e J 
SAeone af 2 ad 6th Mm Oh-tiielsiel eft fo no 
gadis eric Pitt. ©: 4530p ,Sheaditivoes: nolsets lita) = 36, 
1 \Yew vd (wok? twekieve) Mone ener 202 SHctg isoleydg 
i VE ieee SRPOTU Se AOL Tete Gta delevatqged edd pak Cate | 
PAGS 7225080 ,POCeIeSSe noisaesze@ bra «> LeeTqs 
nosy rite: 
“ 7 
x \ é 
A 


a os 


A portion of the supply to the soil profile, normally 
when in excess of its storage capacity, serves as the source 
for positive supply to the Aquifer Storage. This contribution 
is affected by the subsurface physical process of Percolation, 
which is characterized by a saturated flow. In situations where, 
and at times when, gravity flow from the soil profile to the 
water table does not exist, a reverse flow (unsaturated flow) of 
moisture may be effected by Capillarity; the process being initiated 
and maintained by a negative gradient set up by the effect of 
evaporation and/or transpiration at the soil surface. 

The Aquifer Storage serves as a source for the active ground- 
water flow process resulting in Baseflow into the stream. In 
Situations where a local aquifer is in continuum with the expanse 
of a regional aquifer, a portion of the input to the system may 
be accounted for as a loss to deep percolation in the regional 
ground-water flow system. 

Total Runoff, the summation of the Surface Runoff and Baseflow, 
routed by way of channel(s) flow is reproduced at the drainage 
outlet as a Total Runoff Hydrograph. 

Direct losses from the various storages, such as Interception, 
Depression and Detention, Liquid Amount and Snowpack, are affected 
by the process of Evaporation. The Transpiration process by 
plants result in indirect moisture losses from the unsaturated and/ 
or saturated zones. The Evaporation and Transpiration processes 
form the overall Evapo-transpiration process and the sum of these 
losses constitute a negative output from the hydrologic system as 


Actual Evapo-transpiration. 
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Figure 2 of Appendix II, relates to the various instru- 
mentation and/or methods used for quantitative determination 
of the various hydrologic variables or parameters in terms of 
amounts, durations, intensities and frequencies of measurements. 
For example, the input variable, Precipitation, is sampled 
for amount, duration, intensity and areal distribution. The 
instruments used for the sampling being storage gauges and 
recording and standard rain gauges. 

The scheme for each segment in Figure 2 of Appendix II, 
outlines, in general, the standard instruments used for 
measurement, the units of measurement and the frequencies of 
measurement. In cases where there are no standard instruments 
for use in determining the variable or parameter, the most 
widely used methods for carrying out these determinations are 
indicated. For example, surface runoff may be estimated from 
a total runoff hydrograph by applying one of several hydrograph 


separation techniques. 


Component Models 


Figure 3 of Appendix III shows a schematization for component 
modelling, relating briefly some of the types of models which 
may be investigated under each component. A 'model', in this 
context, is used very broadly; that is, it incorporates all 
operations that result in a historical or synthetic data 
sequence of the variable. The text on component models in 


Appendix III consists of a caption for each model and expands 
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the summarized scheme of Figure 3,quoting examples of references 
for each type of model, typical data requirements and brief 
comments on present and projected applications relative to the 
River Basin Research Section. 

The block (Figure 3, Appendix III) for 'Data Management' 
outlines briefly the necessary and prerequisite operations to 
be done on the data prior to the development of a component model. 

A group of component models may be combined, through a 
system of transformation ('Transformation System'), to develop 
an Overall comprehensive catchment model. The generated or 
historical data sequence for each of the component models may 
form part of the data pool for use in tests for the sensitivity 


and performance of the overall model developed. 


An Integrated System Mo del 


Figure 4 of Appendix IV, shows a schematization for an 
integrated (comprehensive) hydrologic system model. Outlined 
are a series of operations through various schemes involving 
the Precipitation input as rain, or excess snowmelt, which is 
routed through the various catchment storages, namely, 
Interception Storage, IS, Surface Storage, SS, Aeration Zone 
Storage, AZS, Interflow (intermediate) Zone Storage, IZS, and 
the Saturation Zone Storage,SZS. Routed out of three of the 
storages are outflow components, namely, Overland Flow, OLF, 
Interflow, INFL, and Baseflow, BFL. These are subsequently 
routed into Streamflow through a linear or non-linear reservoir 


and channel storage or flow translation procedure to produce an 
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outflow hydrograph Total Runoff, TRO, at the basin outlet. 

Some of the basic mathematical expressions that relate the 
various parameters within a given storage function are cited 

for the intermediary schemes of operations between the input 

and output. Several explicit relationships have been developed 
by several investigators and were cited by them to be operational 
in an integrated system model through empiricisms and numerical 
approximations. It should be noted, however, that the parameters 
defined in the various functions have to be determined and 
optimized on the basis of sufficient physical justifications. 

The attached text to Figure 4 of Appendix IV, serves to 
define the various parameters of the given functions and outlines 
further a brief account of the operational sequences throughout 
the system. It is not the intention at this stage to elaborate 
on the utilization of these functional relationships, in their 


explicit forms, in actual numerical or computational procedures. 


A Specific Model Approach 


Figures 5a to 5c of Appendix V, are representations of 
basic schemes for a more detailed and specific approach towards 
a framework for a type of parametric model. Starting from a 
simple concept and going through some basic definitions, the 
model is built up to an operating format. Figure 5a of 
Appendix V, is a definition scheme for a parametric model 
(Diskin, 1970) for use in an interpretation of a linear regres- 


sion relationship between precipitation and runoff. Figure 5b 
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of Appendix V, shows an expansion of the model, taking into 
consideration implicitly, retention storage and evapc-transpiration. 
Figure 5c of Appendix V, is the flow chart for the parametric, 
linear regression model outlined in Figure 5a and 5b of Appendix 

V. Given is the flow of mathematical operations and computations 
that may go into the simulation, assuming that the necessary 
parameters, constants and coefficients (regressions and implied 
basin characteristics) were previously determined from other. 
regression sub-routines or by other estimating methods. The 

text on the lumped-parametric model of Appendix V, outlines further, 
the step-by-step development of the model, with definitions for 

the various parameters, coefficients and constants, with an 
explanation of their implied physical significance. Further 
treatment in the text presents a section on the relative stat- 
istical criteria to be used for simplification of the model (e. g. 
operation in one or two elements) for an appropriate application 


in the regression sub-routines. 


A Brief on Stochastic Models 

Figure 6 of Appendix VI, shows a framework for a stochastic 
hydrologic systems model (Chow and Kareliotis, 1970), outlining 
schematically the inter-relationship (formulated on the principle 
of conservation of mass) among the individual models, each of 
which is a component stochastic model. A very long sequence 
of historical record is required for each time series in order 


to adequately determine its probability distribution, or to 
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determine by correlogram or spectral analysis, the type of 
stochastic model (autocorrelation, harmonic, moving average, 
etc.) to be used in the generating processes. The text of 
Appendix VI gives a short review on an approach to stochastic 


hydrologic system modelling. 
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APPLICATION 

The review on conceptual hydrologic modelling outlines 
broadly the extent to which the River Basin Research Section 
plans to proceed with the formulation (building and calibration) 
of a hydrologic model for each of the Ministry of the Environment - 
IHD representative drainage basins. 

The literature reviewed indicates that not many successful 
attempts have been made to formulate one overall hydrologic 
system model which has a general application to all types of 
drainage basins and hydrologic situations. However, the principles 
and approaches that have been developed from the application of 
modelling to more specific conditions, provide for the selection 
of an approach that may be transposed to areas with similar, basic 
hydrologic characteristics. 

The emphasis in modelling has been on separate hydrologic 
phases (component models). Component modelling, theoretically, 
requires a substantial amount of information and lengthy sequences 
of data for a meaningful approach towards formulating a simulator. 
Also, long sequence of record is required for a totally effective 
calibration of the model. The problem of lack of sufficient data 
and length of records, however, is often alleviated by making 
implicit or explicit assumptions regarding a parameter, group 
of parameters and/or some aspects of the inter-reacting hydrologic 
processes. These assumptions are often desirable for simpli- 


fication of the model to make it computationally tractable. 
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The present extent of the instrumentation in the River 
Basin Research Section's drainage basins and the degree of 
data collection activity, relative to the length of historical 
record of the various hydrologic and geologic variables compiled 
to date, will necessarily dictate the nature and types of 
hydrologic modelling that can be attempted. The relatively 
short sequences of data collected subject the choice of models 
to the confines of parametric and deterministic types and of 
a water balance analysis or accounting approach. Attempts can 
be made as first approximations, at model calibration through 
parameter optimizations, based on averages, ratio estimates or 
least-square methods. 

For the present time, component modelling will be restricted 
to manipulations with the parametric and deterministic types, on 
the basis of tried approaches and existing modelling techniques. 
Knowledge gained from the insight into and/or trials with 
existing model techniques will be used in the eventual formulation 
of an overall comprehensive catchment model, based on the data 
from one or more of the representative basins. 

Although cursory and trial attempts will be made with a 
modelling exercise based on conventional probabilistic treatments, 
no extensive attempts will be made at this stage with stochastic 
modelling due to the lack of extensive data for each component 
time series. Although long series of records are a prerequisite 


for effective utilization of this approach to modelling, attempts 
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may be made to formulate and test simpler existing models, e.g., 
the formulation of a precipitation model generator (Grace & 
Pagleson vao0d/ ss Perkins,.1971)., by utilizing relatively short 
sequences of recording precipitation data to develop a stochastic 
model for generating short-time-interval sequences of rainfall 
data. 

Current activities in the Section are centred around data 
Management (data reduction, compilation,error analysis, etc.) 
which is a prerequisite for a rational and concerted approach 
towards overall water balance analyses and hydrologic modelling. 
Models will initially be of a parametric type, based largely on 
the framework set out in the definition scheme for a parametric, 
linear regression model. Subsequently a more detailed approach 
will be followed, as outlined for an integrated hydrologic 


system model, as shown in Appendix IV. 
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IHD PROJECT REPORT 


Project Title: Representative Basin Studies - East and Middle 
Oakville Creeks. 

File Reference: IWB - RB - 16; ONT - 4 

Objectives: Collection of hydrologic and geohydrologic data 

necessary for the study of all aspects of the hydrologic regime 

and the general water balance in a basin representative of clay 

and clay till plains in southern Ontario. 

Instrumentation: 5 recording and 3 non-recording streamflow 

gauging stations; 1 main meteorological station, 6 satellite 

meteorological stations; 2 recording, ground-water observation 

wells; 21 peizometers in 8 nests, 12 private, non-recording, 

shallow, ground-water observation wells, 9 snow courses and 12 


soil moisture index stations. 
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IHD PROJECT REPORT 


Project Title: Representative Basin Studies - Wilton Creek 

File Reference: IWB —- RB —- 20; ONT -4 

Objectives: poll sation of hydrologic and geohydrologic data 
necessary for investigations into all phases of the general 
water balance in a basin representative of limestone ponerse 
in southern Ontario. 

Instrumentation: 2 recording and 2 non-recording streamflow 
gauging stations; 2 main meteorological stations and 5 satellite 
meteorological stations; 4 recording, ground-water observation 


wells and 11 piezometers in 4 nests. 


iy 

nm = ‘Wh ad 

- ; @ ‘ ar 
68.5 2LVOlO30N 8 oO | foe 

omnt snoléspisagyal a 
q 


7 
_ 
ed 
a 
7 
y« 
: 
7 
7 
7 
= | 


IHD PROJECT REPORT 


Project Title: Representative Basin Studies - Bowmanville, 
Soper and Wilmot Creeks 
File Reference: IWB - RB - 19; ONT - 4 
Objectives: Collection of hydrologic and geohydrologic data 
necessary for the study of all aspects of the hydrologic 
regime and general water balance in a basin representative of 
moraine, till and clay areas in southern Ontario. 
Instrumentation: 16 recording and 2 non-recording streamflow 
gauging stations; 2 main meteorological stations and 12 satellite 
meteorological stations; 12 recording ground-water observation 
wells; 17 piezometers in 8 nests; 23 private, shallow, non- 
recording ground-water observation wells; 16 soil moisture 


index stations; 2 ground-temperature recorders and 12 snow 


courses. 
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IHD PROJECT REPORT 


Project Title: Representative Basin Studies - Venison Creek 
File Reference: IWB - RB - 17; ONT - 4 

Objectives: Collection of hydrologic and geohydrologic data 
HeCecseary 2O0r the study of ‘all aspects of the hydrologic regime 
and general water balance in a basin representative of sand 
plain conditions in southern Ontario. 

Instrumentation: 3 recording and 1 non-recording streamflow 
gauging stations; 1 main meteorological station and 3 satellite 
meteorological stations; 4 recording, ground-water observation 


wells and 3 piezometers in 1 nest. 
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IHD PROJECT REPORT 


Project Title: Representative Basin Studies - Blue Springs Creek 
File Reference: IWB - RB - 18; ONT - 4 

Objectives: Collection of hydrologic and geohydrologic data 
necessary for the study of all aspects of the hydrologic regime 
and general water balance in a basin representative of kame, 

till and bedrock conditions in southern Ontario. 

Instrumentation: 3 recording and 6 non-recording streamflow 
gauging stations; 5 recording ground-water observation wells; 

11 piezometers in 3 nests and 8 private, non-recording, shallow 


ground-water observation wells. 
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HYDROLOGIC MODEL: II-1 
Basic Conceptualization II-1 
Instrumentation II-2 
Figure 1: A Schematization for Functional 
and Inter-related Hydrologic 
Processes II-4 


Figure 2: 


A Schematization for the 
Superposition of the Basin 
Instrumentation Upon the 
Functional and Inter-related 
Hydrologic Processes II-5 
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HYDROLOGIC MODEL: 


Basic Conceptualization 


Figure 1 shows the schematization outlining the functional 
and inter-related processes based on the conventional concepts 
of the hydrologic cycle. 


The end product of the meteorological process, Precipitation, 
1s rain and/or snow, a portion of which is intercepted by the 
vegetal canopy as Interception Storage and eventually lost to 
Evaporation (or Sublimation). Precipitation as a liquid is held 
momentarily at the ground surface as Detention Storage. Some, 
through the process of Infiltration, enters the Unsaturated Zone 
and that in excess of the soil moisture capacity goes to replenish 
the local ground-water flow system (Saturated Zone) by Percolation. 
Precipitation with rates in excess of the infiltration capacity 

at the soil surface contributes to Total Runoff by way of direct 
Surface Runoff. A portion of the precipitation to direct Surface 
Runoff is trapped as Depression Storage in ponds, swamps, etc. 
Subsequently, this retained portion of the runoff contributes 
respectively to Evaporation, to Infiltration to the Soil Profile 
and recharge to the local Aquifer. Portions of the recharge 
eventually contribute as Baseflow to Total Runoff and to Regional 
Flow as Deep Percolation. Eventually, there may be a reverse 

flow to the Unsaturated Zone as a result of a capillary (suction) 
potential set up by Evaporation and/or Transpiration at the soil 
surface. 


Precipitation which occurs as snowfall is stored on the 
ground as Snowpack. The pack undergoes various forms of met- 
amorphosis prior to the active snow-melting process. During 
the metamorphosis phase, a quantity of moisture is lost to 
Evaporation (by Sublimation); also during this period, several 
phases of heat exchange take place at the snowpack-soil interface 
and between the snowpack and the liquid rain. These and several 
other types of energy fluxes (Radiation, Temperature) influence 
the onset of the Snowmelt process by way of convective heat 
transfer between the pack and atmosphere. The Snowmelt, as 
liquid excess drained from the Snowpack, follows a similar path 
to that of the rain; that is, via Depression and Detention 
Storages and Infiltration, and subsequently to Surface Runoff 
and Baseflow to Total Runoff which are reproduced as a Total 
Runoff Hydrograph at the basin outlet. 


Transpiration by crops from the Unsaturated and/or 
Saturated zones, together with evaporation from the various 
processes (Interception, Snowpack, Detention and Depression 
storages), contributes to the Actual Evapotranspiration, 


from the basin. 
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The ancilliary meteorological variable Temperature, T, 
is associated with Radiation, R, (as longwave, shortwave or net 
radiation), with net energy (heat) flux at the incident surface. 
The Humidity, H, is a function of the vapour pressure, a and 
the saturated vapour pressure, es, at a given air temperature. 
The Wind, W, is expressed in terms of the wind speed, u, and 
direction. Sunshine duration, S, is directly related to the 
temperature and shortwave radiation variables. These ancilliary 
meteorological variables are necessary prerequisites for use in 
computation (through empirical formulae) of snowmelt and eva- 
poration from snowpack and evaporation from available sources of 
free water. 


Instrumentation 


Figure 2 shows the various instrumentations and/or methods 
used for quantitative determination of the various hydrologic 
variables or parameters. 


The rainfall is measured by standard rain gauges, RS, 
providing total depth quantities in inches, arranged to yield 
daily, monthly and yearly summaries. Rainfall intensities 
and durations are measured by recording rain gauges, RA, providing 
summaries of storm amount, duration and rates per unit time which 
are the data necessary for such analytical procedures as storm 
intensity-frequency and depth-area-duration analyses. Snowfall 
may be measured accumulatively with Alter-shielded or Nipher- 
shielded storage gauges; accumulated snowfall on the ground as 
snowpack is core-sampled for depth, water equivalent and density. 


The other meteorological variables are measured by con- 
ventional instruments: temperature extremes by maximum and 
minimum thermometers; humidity and temperature by hygrothermo- 
graphs; duration of sunshine by Campbell-Stokes solarimeters; 
longwave and shortwave (Solar) radiation by Eppley or Kipp 
pyranometers; net radiation by CSIRO pyrradiometers; evaporation 
by Class-A pans and wind speed by anemometers. 


Interception (amount of precipitation intercepted by the 
vegetation) may be determined by plot studies, or estimated with 
the use of existing, general empirical equations. 


Snowmelts are usually determined by computations from 
snowmelt indices and/or with general empirical equations. In 
the synthesis of snowmelt hydrographs, some basic factors have to 
be determined. These are areal extent of snow cover; snowpack 
parameters with respect to depth, water equivalent, temperature 
and free water; the variation of these parameters with elevation; 
albedo of the snow surface; and the critical sequence of factors 
affecting melt such as melt rate and loss and runoff conditions 


of the basin. 
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The extent of the depression storage capacity and its 
distribution may be determined from topographic maps and field 
surveys of the drainage area. 


The infiltration function may be approximated by empirical 
procedures such as sample plot studies (infiltrometer, rain- 
ulators), or by laboratory experiments with field soil samples. 


A time series of soil moisture data may be determined by 
thermogravimetric analyses of field samples or monitored in 
Situ by tensiometers or electrical resistance block methods or 
by the nuclear (neutron meter) method, etc. Soil temperature 
profiles may be monitored in situ by thermohms or by ground- 
temperature recorders. The necessary soil moisture characteristics, 
(soil moisture-hydraulic conductivity and soil moisture-tension 
curves), for describing the unsaturated flow, may be determined 
through laboratory experiments on specific soil samples. 


Aquifer characteristics (hydraulic conductivity, transmis- 
Sibility) may be determined by pumping tests and by analytical 
techniques using ground-water level fluctuations in observation 
wells and piezometers. The saturated hydraulic conductivity 
characterizes the rate of percolation. 


Total runoff may be determined by the use of conventional 
streamflow gauging stations (recording gauges with natural 
controls, NR, or artificial controls, AR and staff gauges, NN). 
Measured stages at each gauge are converted to continuous 
discharge as mean daily flows in cubic feet per second, with 
the use of empirically developed stage-discharge rating curves. 


Baseflow and direct surface runoff may be determined by 
analytical methods such as hydrograph separation. A direct 
estimate of baseflow from some areas may be determined by the 
measurement of spring discharges. 


Estimates of the potential evapotranspiration may be 
obtained from adjusted evaporation pan measurements or estimated 
with the use of existing empirical formulae. Estimates of actual 
evapotranspiration may be obtained from lysimeter studies or 
from adjusted potential evapotranspiration with a crop consumptive- 
use factor determined from land-use studies. 
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Precipitation 


storage gauges, standard 
rain gauges, recording 


gauges, etc. 

inches, ins /hour, ins /day 
continuous, hourly, daily, 
year around 


Interception 


- rational estimates 
(empirical equations 
based on plot studies): 


Rain 

I: standard & recording rain 
gauges, RS, RA 

UM: inches, ins /hour, ins /day 

F: continuous, hourly, daily, 

year around 


I: 


F: 


Zecorders 
UMS, 3G 


Snowfall 
storage gauges,snow board 
UM: inches, ins /day 

intermittent - snow season 


Conductive Heat Exchange 
I: thermohms, ground- temperature 


F: intemittent - snow season 


Sunshine 


I: solarimeter 
: hours/day 
continuous, daily, 
year around 


I: 
UM: 
F: 


I: 


Infiltration 
infiltrometer, rainulator 
ins /hr. 

sample plot studies 


Soil Profile - Unsaturated Zone 
thermogravimetric, electrical resistance 

blocks, neutron meters, network or plot studies 
ins, % moisture ; 
intermittent - non-winter months (daily, 
weekly or monthly surveys) 


Temperature 


I: maximum-minimum 
thermometers 

um: °F, °c 

F: daily, year around 


Snowpack 

I: snow course, snow stakes, 
snow pillow 

UM: depth, water equivalent, 
density, temp., areal 
distribution 

F: intermittent - snow season 

(weekly or bi-weekly surveys) 


pyranometer, 
pyrradiometer 2 
Langley (cal./cm /min.) 
continuous, tntegrated 
hourly, year around 


Snowmel t 

I:] melt plot, snow pil) 
UM: ins /hour, ins /day 
F; intermittent - snow 


Capillarity 
(unsaturated flow) 

unsaturated hydraul: 
conductivity, soil | 
moisture-tension | 
by lab. determinatic 


Recharge ~ Per colati 
(saturated flow) 


saturated hydraulic — 
tivity and total pot 
head 

field studies 


Wind Evaporation 


anemaneter 
UM: mi/hour, mi/day 

F: continuous, integrated hourly, 
year around 


I: Class 'A' pan & auxiliary 
instruments 

UM: ins/day 

F: intermittent - non-winter 

months 


NOTE: I =- Instrument 
UM - Units or Measurements 
F - Frequency of Measurements 


Evapotranspiration 

I: evaporation pan, evaporimeter, 
lysimeter 

UM: inches, ins/day 

F: intermittent - non-winter 

months 


1 
tive use 


Depression Storage 


(pond, swamps, etc.) 
% - areal estimate 


ter-use 


; 


Surface Runoff 
- computed estimate from 
hydrograph separation 


Total Runoff 
I: streanflow gauging 


stations, AR, NR, AN, NN 


cfs, gpm, gpd, ac -in., 
ac -ft., etc. 


continuous, year around 


Aquifer - Saturated Zone Baseflow 


at 


observation well, piezometric 
nests, GR, P (N). 
foot, hydraulic conductivity, 


Computed estimate 
(hydrograph separation) 
- spring discharge 


transmissibility, potential gradient, 
specific capacity, etc. 
continuous, daily, weekly, year around 


Figure 2: HYDROLOGIC MODEL: A Schematization for the Superposition of Basin 
Instrumentation upon the Functional and Inter- 


related Hydrologic Processes 
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APPENDIX III 


HYDROLOGIC MODEL: 
Component Models 


Precipitation Models 
Method 
Data Requirement 
Projected RBR Application 


Snowmelt Models 
Method 
Data Requirement 
Projected RBR Application 


Soil Moisture Models 
Method 
Data Requirement 
Projected RBR Application 


Evapotranspiration Models 
Method 
Data Requirement : 
Projected RBR Application 


Surface Runoff Models 
Method 
Data Requirement 
Projected RBR Application 


Ground-Water (Baseflow) Models 
Method 
Data Requirement 
Projected RBR Application 


Total Runoff Models 
Method 
Data Requirement 
Projected RBR Application 


Transformation System (Catchment) 
Models 

Method 

Data Requirement 

Projected RBR Application 


Figure 3. A Schematization for 
Hydrologic Component Models for 
Data Generation and Simulation 
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HYDROLOGIC MODEL: 


Component Models 


In general, a component model refers to a mathematical 
formulation, involving operations with a data sequence of a 
single hydrologic process. This may be used for extending the 
data sequence or used to produce synthetic quantitative traces 
of the variables for use in simulation. These component models, 
schematized in Figure 3, are treated on the basis of lists of 
different methods and existing procedures that may be used for 
a model development. Quoted are some of the references relating 
to types of models. 


The likely applicability of a type of model(s) to River 
Basin Research (RBR) present and projected investigations in 
the representative basins is also briefly discussed. 


The cited references are given as examples of different 
types of investigations relating directly or indirectly to 
modelling. An attempt was made to identify the reviewed 
literatures and group them under the different component models. 
Identification was based on the type of hydrologic process and 
the choice of the method used in the modelling. 


Precipitation Models 


Method: Deterministic 
i) inter-station relationship (correlation and 
regression models) for use in data supple- 
mentation and extension. 
ii) storm transposition. 


Probabilistic 

i) order statistics, e.g. conventional procedures 
for depth-area-duration analysis. 

ii) frequency analyses (probability of discrete 
events), e.g. storm intensity -frequency 
analysis and extreme values analysis. 

iii) random sequences, (stationary and non- 
stationary process), e.g. stochastic models 
for synthesis of hourly rainfall data, models 
of the types formulated by Grace and Eagleson 
(1967)*, Pattison (1965), Franz (1971), 
Raudkivi and Lawgun (1970), Hiemstra and 
Creese (1970). 


* References appended in Bibliography. 
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Data Requirement 


historical data: daily, monthly and yearly summaries from 
the basin network stations; recording precipitation gauge 
data summaries; monthly, yearly, summaries from selected 
regional index precipitation stations. 


Projected RBR Application 


each basin has 1 or 2 recording rain gauges and a network 
of standard rain gauges which provide the required data 
summaries. 


historical records of several long-term regional stations 
within each climatic region are available from monthly 
climatological station reports, published by the Atmospheric 
Environment Service, Canada. 


application in modelling may be attempted with data extensions, 
order statistics and sample computations based on frequency 
analysis and probability of discrete events. Some attempts 
may be made with stationary stochastic methods. This involves 
utilizing a time series of recording rain gauge data to 
formulate a stochastic model for use in simulating short- 
sequence rainfall data. 


' Snowmelt Models 


Method: Deterministic 

i) empirical equations - snowmelt synthesis 
with the use of general empirical snowmelt 
equations of the Corps of Engineers, U. S. 
Army (1960). 

ii) melt indices - use of melt indices, e.g. 
degree-day index, or melt-plot studies, model 
byeWatt andiHsu (297). 

iii) index plot and empirical equations - 
snowmelt simulation using a hybrid of index 
plot studies and empirical formulae, e. g. 
models formulated by Anderson and Crawford 
(1964), Amorocho and Espildora (1966), 
Crawford and Linsley (1966). 


Data Requirement 


- basic meteorological parameters (rainfall, temperature, 


radiation, vapour pressure, duration of sunshine, wind, 
convective heat transfer). 


areal extent of snow cover and snowpack parameters (depth, 
water equivalent, density, temperature, albedo of snow 
surface, permeability, etc.). 


streamflow records. 
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Projected RBR Application 


- basic meteorological parameters are monitored continuously 
in each basin. 


- areal extent of snow cover and snowpack parameters (depth, 
water equivalent, density, etc.), are determined seasonally 
in two baSins. 


- winter streamflows are monitored continuously. 


- present application in modelling may be attempted with 
snowmelt hydrograph synthesis utilizing general melt 
equations and/or melt indices, e. g. models of the type 
formulated by Watt and Hsu (1970) which utilize a melt- _ 
index based on the degree-day concept, or a modified approach 
to the models used by Amorocho and Espildora (1966) which 
utilize empirical equations based on energy-melt concepts 
and index plot observations. 


Soil Moisture Models 


Method: Deterministic 

i) *lumped-system as related to conceptual 
catchment storage; the use of soil moisture 
data, for specified soil types, in a budget- 
type equation to develop empirical relation- 
ships for use in estimating soil moisture 
regimes, e.g. models of the type used by 
Zahner (1967), Visser (1966), Baier (1969), 
Baier and Robertson (1965). 

- the use of depletion curves and accretion 
regressions from measured field data to 
develop empirical relationships for predicting 
moisture regimes, e. g. models of the types 
used by Carlson et al (1956), Stearn and 
Carlson (1960). 

ii) lumped-system as related to infiltration; the 
use of soil moisture parameters in the 
modelling of infiltration process with methods 
such as those used by Overton (1964), 

Huggins and Monke (1968), Deboer and Johnson 
(1969), Skaggs et al (1969), Ligon, et al (1965). 


Data Requirement 


- soil types, land use and cropping practices and their areal 
distribution. 


- hydrologic characteristics (moisture retention curves, 


* "lumped-system' - the spatial distributions of the input and 
output are ignored. 
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infiltration rates, hydraulic conductivity, etc.) of the 
major soil types. 


- historical data of absolute soil moisture or changes in 
Soll moisture for the study area, point measurements from 
a network of stations. 


- meteorological variables (rainfall, temperature, radiation, 
humidity, evaporation, etc.). 


Projected RBR Application 


-~ soil moisture observations (absolute and changes in moisture 
from field studies ) are being carried out in one sub-basin 
(BSW-Wilmot) and the initiation of a similar study is in 
progress for EMO. 


- summaries of major soil types and land use and cropping 
practice in these basins are compiled. 


- soil hydrologic characteristics to be determined by lab- 
Oratory experiments on field samples. 


- present application to modelling may be attempted with 
deterministic methods, such as the use of a water balance 
type equation to develop empirical relationships from field 
observations for use in estimating soil moisture regimes, 
e.g. method by Zahner (1967), Baier and Robertson (1965) ; 
the use of correlation relationships between soil moisture 
changes and meteorological variables to predict moisture 
regimes, e.g. method by Stearn and Carlson (1960); and the 
use of soil moisture parameters in infiltration models, e.g. 
method of Huggins and Monke (1968). 


Evapotranspiration Models 


Method: Deterministic 
4.3) empirical relationships for potential evapo- 

transpiration; the use of developed evaporation 
formulae which are based on the Energy-Budget 
and Mass-Transfer theories of evaporation such 
as the Penman or Thornthwaite formulae and 
method by Gay (1970); the use of adjusted 
Class-A pan data with other meteorological 
variables to develop correlation or regression 
relationships for use in predicting potential 
evaporation regimes, such as done by Mukammal 
and Bruce (1960), Christiansen (1966), Morton 
(1971), Harmon (1966). 
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ii) empirical relationships for determining actual 
evapotranspiration; the use of a water balance 
type equation with measured field data to 
develop an empirical relationship for use in 
predicting the evapotranspiration regime, 
such as the method used by Bowman and King 
(1965); the use of a meteorological budgeting 
technique with meteorological and soil moisture 
data to develop empirical relationships for 
estimating daily evapotranspiration, e.g. 
method by Baier (1966); the use of pan eva- 
poration and meteorological variables to 
develop relationships for estimating actual 
evapotranspiration such as summarized by 
Hargreaves (1960); the use of solar radiation 
and temperature to estimate evapotranspiration, 
e.g. method by Robb (1966). 


Probabilistic 
1) time-series model; the use of time-series 
analysis on historical sequences of estimated 
potential evapotranspiration to determine a 
stochastic model for data generation and 
Simulation, Chow and Kareliotis (1970). 


Data Requirement 


- historical data: hourly, daily, monthly or yearly summaries 
of meteorological data (radiation, temperature, humidity, 
wind speed, duration of sunshine, Class-A pan data, etc.). 


Projected RBR Application 


- the basic meteorological variables including Class-A pan 
evaporation data are monitored in each basin. 


- application to modelling may be attempted with computations 
based on available empirical formulae or by trial computation 
with formulations through correlation and regression analysis, 
e. g. the application of a modified Penman's or Thornthwaite's 
equation and the development of regression relationships 
between Class-A pan evaporation and meteorological variables 
for use in data supplementing or data extension, or the use 
of a budget equation with measured meteorological and soil 
moisture data to develop relationships for use in estimating 
actual evapotranspiration. 


Surface Runoff Models 


Method: Deterministic 
i) linear time-invariant; the use of a distributed 
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linear method for the synthesis of runoff 
hydrographs from given rainfall inputs, such 
as the Unit Hydrograph Method, as applied by 
Eagleson (1962), Mitchell (1967), Gray (1961), 
Viessman (1968), Overton (1967), Dooge (1959); 
Eagleson et al (1965); or the use of other 
routing methods such as the rational method 
and the summing of hydrographs from a 

number of sub-areas within a basin. 

~ the use of the physiographic parameters of 
the basin and surface runoff characteristics 
to develop regression equations for use in 
Simulation, e. g. as developed by Narayana 
and Riley (1971), Schulzgvet lai (197i), 

Taylor and Schwarz (1952), Larson (1965). 

- the use of a hydraulic model of the catchment 
to simulate runoff, e.g. based on Kenematic 
Faves cengty as applied by Wooding (1965, 
1966). 

ii) non-linear time-invariant; the use of a 
distributed non-linear method to synthetise 
runoff hydrographs, e.g. IUH as applied by 
Singh (1964). 


Probabilistic 
i) frequency analyses; the use of peak flow 
data for annual or partial duration series 
to develop frequency curves and statistics 
for use in flood flow predictions, e. g. 
as applied by Snyder (1958). 

- the use of the statistical parameters of 
flood series and rainfall extremes to develop 
regression relationships for use in simulation, 
e. g. aS applied by Reich (1970), Shane and 
Gaver (1970). 

ii) random sequences (stationary and non- 
stationary processes); the use of the 
stochastic property of flood flow series to 
develop a model for use in simulation, e.g. 
as applied by Todorovic and Zelenhasic 
(1970), Shane and Lynn (1964); the application 
of a stochastic model to storm runoff and 
rainfall data to formulate relationships 
for use in generating sequential data, 
after Chow and Ramasseshan (1965). 


Data Requirement 


- historical data: streamflow (mean daily and peak flows); 
recording precipitation summaries (amount, intensity, duration). 
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Projected RBR Application 


- the use of deterministic and probabilistic methods - re 
Synthesis of surface runoff hydrograph, with the use of unit 
hydrograph or instantaneous unit hydrograph techniques; re 
development of relationships between runoff hydrograph 
parameters and basin physiographic parameters; re development 
of frequency curves and parameter relationships for flood 
flow series; re development of regression relationships between 
flood flow and rainfall extreme series, etc. 


Ground-Water (Baseflow) Models 


Method: Deterministic 


i) 


lumped-system; to consider the relationship 
between total quantity of precipitation and 
its direct runoff component, for use in 
estimating ground-water flow and rainfall 

loss relative to the infiltration capacity 
and soil moisture as applied by Kadoya (1967), 
Meyboom (1961). 


ii) *distributed-system; to consider the theory of 


flow problems re models based on the analytical 
and numerical solution of the flow equation with 
given boundary conditions for specified problems, 
e. g. for a simplified, two-dimensional flow 
problem as applied by Toth (1962, 1963) and 

for multi-layered, two-dimensional flow 

problem as applied by Freeze and Witherspoon 
(1967); to consider heterogeneity of basin 
models by determining the statistical variability 
of input data (conductivity, specific yield, 
water-level elevation, spatial dimension of 
aquifer, etc.) e. g. solution of the flow 
equation with variable coefficients as applied 
by McMillan (1966), Shahbazi and Todd (1967), 


Kivaly? (19 7714)2: 


Probabilistic 


i) 


correlation; to develop statistical relation- 
ships between baseflow data (determined by 
hydrograph separation techniques) and drainage 
basin parameters or meteorological parameters, 
e. g. relationship between drainage density 

and baseflow as applied by Trainer (1969), 
between baseflow and basin parameters as applied 
by Furness and Busby (1969), between ground- 
water chemistry and baseflow as applied by 


* 'distributed-system' - the spatial distribution of the input 


and output is considered. 
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Pinder and Jones (1969); relationships of 
quantitative geomorphology to baseflow as 
applied by Farvolden (1963); relationship 

of recharge to ground water from storm inputs 
to increment in baseflow as applied by 
Trainer (1969). 


Data Requirement 


- surficial and subsurface geology. 

- summaries of water-table fluctuations, spatial dimension of 
poe chemical properties of ground water, pumping tests, 
etc. 

- aquifer characteristics and their variability (specific yield, 
hydraulic conductivity, etc.) 

- natural spring discharge measurements (where possible). 

- low-flow data (streamflow hydrograph) . 


Projected RBR Application 


- compilation of basic data (as above) is in progress. 

- the use of parametric methods such as the lumped-system to 
relate indices of gross precipitation inputs, direct runoff 
components and corresponding antecedent moisture conditions. 

- the use of simplified one or ii aumpapeendecis 2s 32s flow problem 
models. 

- the use of correlation relat tonani os between baseflow components 
and drainage basin or meteorological parameters. 


Total Runoff Models 


Method: Deterministic 
in) linear time-invariant; the use of an 
antecedent-precipitation-index type of 
rainfall-runoff relationship and a 
retention index to simulate continuous 
streamflow hydrographs from rainfall inputs, 
e.g. as applied by Sittner et al (1969). 


Probabilistic 

ds) regression and correlation; the use of 
multiple regression techniques to develop 
relationships between streamflow data at 
inter-related stations; to develop relation- 
ships between streamflow parameters and 
basin physiographic parameters; to develop 
relationships between the original variables 
or statistical parameters of rainfall and 
streamflow series, as applied by Beard (1965), 
Benson and Matalas (1967), Caffey (1965). 
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ii) auto-correlation; the use of stochastic 
models (Markov Models) on historical stream- 
flow data to develop relationships for use in 
generating synthetic data; the use of a 2nd- 
Order Markov model, applied to standardized 
variates (assumed log-Pearson Type III 
distribution), e. g. method by Beard (1967); 
the use of Markov model for low-flow analysis, 
e. g. method by Fiering (1964); the use of a 
lst-order Markov model and regression analysis 
on multiple-station flow data to develop a 
Simulator for generating daily flows, after 
Payne and Newman (1969). 


iii) other time-series analyses; the use of time 
series analyses on historical streamflow data 
to determine the inherent stochastic model for 
use in data generating; the autocorrelation 
and spectral analysis of streamflow series, 
after Quimpo (1967), Roesner and Yevdjevich 
(1966), Svandize (1967). 


Data Requirement 


- historical data: streamflow (mean daily, monthly and yearly 
summaries) for both base and secondary gauging stations; 
precipitation summaries. 


Projected RBR Application 


- compilation of streamflow data for each basin gauging network 
is in progress. 

- application of deterministic and probabilistic methods - re 
synthesis of total runoff hydrographs, using a developed 
relationship between antecedent-precipitation-index, retention- 
storage-index and runoff, or using a multi-station regression 
relationship, or using a regression between streamflow and 
basin parameters and/or meteorological parameters. 

- application of probabilistic methods - re use of simpler Markov 
models for data generation (generating models based on mean, 
standard deviation and serial correlation of a flow sequence). 


Transformation System (Catchment) Models 


Method: Deterministic 
me) lumped-system; to consider the catchment as a 


linear system, governed by a differential 
equation, the integration of which provides 
the convolution equation for the linear 
transformation of an input (index of precip- 
itation) to an output, as applied by Amorocho 
and Orlob (1961), O'Connell and Nash (1970). 
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distributed-system; to consider the catchment 
to be composed of an infinite array of lumped, 
independent. linear subsystems with each sub- 
system being governed by a linear differential 
equation, the summation and integration of 
which leads to the convolution relationship 
for transforming inputs to outputs, as applied 
by Eagleson (1967), Wooding (1965, 1966), 
Brakensiek and Onstad (1968), Huggins and 
Monke (1970), Foster, et al (1968). 


linear reservoir; to consider the catchment as 
a Single reservoir or a cascade of n-linear 
reservoirs, in each of which the storage is 
linearly related to the outflow, e. g. use of the 
storage equation and the hydrologic continuity 
equation to derive the differential equation 
that governs the linear reservoir, as applied 
by Nash (1957), Sarma et al (1969), Overton 
(1967)>FOtConnelis etl aleil970). 


linear channel and linear reservoir; to 
consider the catchment to consist of a series 
of linear channels and linear reservoirs 
having a combined linear translation effect, 
using a time-area-concentration diagram as 
applied by Dooge (1967). 


non-linear channel and reservoir; to consider 
the catchment storage to be distributed and 
non-linear by dividing the catchment into 
Sub-areas, considering their individual 
inputs and routing the excess from each sub- 
area through different amounts of storage, as 
applied by Laurenson (1964), Deboer and Johnson 
(1969); to consider the operation of the 
catchment by a functional series (the gen- 
eralization of the U.H. concept to that of a 
non-linear response function) as treated by 
Amorocho and Orlob (1961), Machmeier and 
Larson (1967), Overton (1967), Prasad (1971). 


Probabilistic 


i) 


di) 


the use of a linear regression relationship 
between inputs and outputs of the system, as 
applied by Diskin (1970). 


to consider the catchment as a stochastic hydro- 
logic system; formulate a system model on the 
basis of conservation of mass regarding the 
components (precipitation, runoff, storage 

and evapotranspiration) as individual stochastic 
processes, aS applied by Chow and Kareliotis 
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Data Requirement 


physiographic parameters of the drainage basin (area, 

drainage density, slope, main stream length, etc.). 

Ree data or hydrologic parameters (precipitation, runoff, 
etc. 


RBR Projected Application 


the use of the lumped-systems re application of a convolution 
integral to transform lumped inputs to outputs. 

the use of simpler distributed system re application of the 
convolution integral to sub-systems. 

the use of a single or double linear reservoir storage re 
application of the storage and hydrologic continuity equations 
to derive the differential equation that govers the input- 
output relationship. 

the use of linear channel translation and linear reservoir 
storage re application of a time-area-concentration relationship 
for the transformation of inputs to output. 

the use of a linear regression relationship between precipitation 
and runoff. 
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HYDROLOGIC MODEL: 


Figure 3. A Schematization for Hydrologic 
Component Models for Data 
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APPENDIX IV 


HYDROLOGIC MODEL: 


An Integrated Hydrologic System Model 
Precipitation 
Interception Storage 
Surface Storage 

Depression Storage 

Detention Storage 

Infiltration 
Aeration Zone Storage 
Interflow Zone Storage 
Saturated Zone Storage 
Total Runoff 
Evapotranspiration 
Working Model 


Figure 4: A Schematization for an 
Integrated (Comprehensive) Hydrologic 
System Model 
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HYDROLOGIC MODEL: 
An Integrated Hydrologic System Model 


Figure 4 shows a schematization for an integrated hydrologic 
system model. 


Precipitation 


_, Precipitation, P, aS rain and/or snow (snowmelt) is the 
initializing input to the system. The initial amount is inter- 
cepted to some degree by the vegetal canopy. 


Interception Storage 


The intercepted precipitation is held by the vegetal canopy 
as Interception Storage, IS. This Bis a function of the amount 
of P, the duration of the precipitation, ths the evaporation 


rate, E! at the intercepting surface, the interception storage 


capacity, ISC, of the intercepting area and the ratio of the 
evaporating surfaces to the total horizonally projected canopy 
Surface area, A. An explicit function relating these parameters 
is required for use in determining IS from the gross P. 


Surface Storage 


The excess precipitation goes to satisfy the Surface Storage, 
SS, requirement. The Surface Storage is composed of three 
conjunctional functions, namely the Depression Storage, DS, the 
Detention Storage, D, and Infiltration, I. Some explicit function 
is required to distribute the excess precipitation to the require- 
mentsof DS, D and I. Generally I reacts simultaneously with DS 
and D to the precipitation input, but DS is normally satisfied or 
partially satisfied before the initiation of D. 


Depression Storage 


The Depression Storage, DS, is a function of the depression 
storage capacity, DSC, of the drainage area, the amount of 
precipitation excess, Pays the infiltration characteristics of the 


soil, ic, and the ratio of the depressional storage area to the 
total drainage area, B. 


Detention Storage 


The Detention Storage, D, is a momentary storage which is a 
function of the overland flow length, L, the infiltration 
characteristics, ic, and the overland flow rate, q. The flow 
rate q, in turn, is a‘function of the overland flow length, L, 
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the slope, S, and the flow characteristics (Manning's roughness 
coefficient, n) of the overland flow path and the rate of 
precipitation input, i. Whenever the input rate, i, is greater 
than the rate of infiltration, overland flow, OLF, is initiated 
and is routed out of Detention Storage D in terms of as 


The OLF out from D is routed to the streamflow according 
to some method of linear or non-linear channel storage and flow 
translation to produce the direct Surface Runoff, SRO, component 
of the Total Runoff, TRO, measured at the basin outlet. 


Infiltration 


Tne@inisitration rate, I, is a function of the \soil, mosture, 
SM, the amount of excess precipitation, Par the rate of precip- 


Weattonmtmouc,, 1, the initial infiltration rate, Toe and the 
infiltration capacity, Tye of the soil. Infiltration proceeds 


at a decreasing rate with time; the initial rate is dependent 
on the antecedent soil moisture condition and continues at a 
rate subject to Tue until the liquid amounts in DS and D become 


limiting. Several empirical infiltration formulae have been deve- 
loped by many investigators (Horton, Holtan, Huggins and Monke, 
Phillip, et al). A few of these empirical formulae have been 
Suitably applied to describe the infiltration process in source 
accounting type hydrologic models. The infiltration rate function 
may be examined by the standard flow equation which relates the 
inflow rate, NOS to the ‘hydraulic conductivity, ‘kK, and the ditf- 


fusivity, (D = K/ d (SM),), of the soil and to the differential 
dy 
change in soil moisture with depth, 9 (SM),relative to the soil 
5 P) 


surface. An empirical determination of K and D and a numerical 
approximation technique for solution, are required for incor- 
poration of the infiltration rate function into a working model. 


Aeration Zone Storage 


Infiltrated water moves through the soil profile in 
accordance with known physics of flow through porous media 
thereby satisfying the moisture requirement of the Aeration 
zone Storage, AZS. The soil moisture, SM, is a function of the 
flow characteristics of the soil, i.e. unsaturated hydraulic 


COUGUCTIVA LY, A ais at , capillary potential,w, gravity potential, 


Viramad total porosity, ae of the medium. Flow of water in the 


unsaturated zone can be described by the expanded partial 
differential Darcy flow equation, which is the same standard 
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flow equation for infiltration above, but written in the more 
classical form, 9 ML) = Kngate + 1), relating the rate 


of change of moisture content with the hydraulic potential. 


Interflow Zone Storage 


Infiltrated water, drained by gravity from the Aeration 
Zone Storage, AZS, percolates to the Interflow Zone Storage, IZS. 
Saturated flow is characteristic of the IZS zone. This zone is 
in fact a transient saturated zone or a saturated zone created 
by a less permeable layer below the more permeable soil profile, 
thereby forming a 'perched water table'. The physics of 
Saturated flow can be used to describe the seepage to be routed 
out of this IZS, as interflow, INFL, which is subsequently routed 
into the streamflow as part contribution to Total Runoff, TRO. 


Saturated Zone Storage 


Drainage by percolation from the IZS goes to recharge the 
Saturated Zone Storage, SZS. Moisture flow in this zone is 
characteristically saturated flow. The outflow q, is dependent 
on the amount Q5 in storage and the total hydraulic potential, 


H; the outflow rate being regulated by the saturated hydraulic 


conductivity, rane For saturated flow, SM and A atic are 


constants; and 9(SM)/dst = o. With these known conditions the 
differential flow equation may be solved for unconfined steady- 
state flow in an isotropic porous medium resultingwith the Laplace's 
equation V(KVH) =O. For a given drainage problem with defined 
boundary conditions, the differential flow equation for the 
problem may be expressed explicitly. With the use of numerical 
approximation techniques, the explicit flow equation may be 

solved approximately and so facilitates the routing of a base- 
flow component BFL, out of SZS, into streamflow as another part 


Poneri1oucion. to TRO. 


The sum of INFL and BFL constitute the subsurface runoff, 
SSRO, component to TRO. 


Total Runoff 


The Total Runoff, TRO, may be routed through the channel 
reach by a linear (or non-linear) storage equation utilizing 
the law of continuity in the basic routing procedure. This TRO 
may be reproduced at the basin outlet as the streamflow hydro- 


graph. 
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Evapotranspiration 


Actual Evapotranspiration, EY, 1s directly proportional 
to the amount of moisture available in the various storages, 
SSM, to the amount of available energy fluxes in terms of 


radiation, R, and temperature, T, to a crop use factor Ce and 


wind speed u, and is indirectly proportional to the relative 
humidity, H, of the environment. 


Evaporation takes place, from the IS and SS at the potential 
rate until the free water available in the sources is used Up . 
Evaporation also takes place from AZS, IZS and/or SZS, ata rate 
which is proportional to the rate at which the water is supplied 
tothe. soil .surface through capillarity. 


Transpiration, by plants, from the AZS, IZS and/or SZS is 
controlled by the amount of moisture available to the plant. 
im the casewor the AZS, this availability of moisture is Jimited 
tO the (range between that of field capacity ‘and wilting’ point. 


An explicit function (not shown diagrammatically) such as 
One relating ET to the potential evapotranspiration rate and soil 
moisture, SM, in the AZS, to crop-use, is required for incor- 
poration into the model to account for ET loss from the basin. 


Wo rking Model 


The development of a working model based on, or modified 
from, Figure 4 requires that a choice be made of an existing 
explicit function for each component of the model. Explicit 
functional relationships may also be developed independently 
through empiricism based on observations and plot studies. A 
linkage or step relationship is also required for the computational 
sequences between the different component schemes. 
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LIST OF SYMBOLS IN MODEL 


IS = INTERCEPTION STORAGE - in. 


P = gross precipitation - in. 

tL. oo duration of stom - in/hr. 

ES = rate of evaporation from the intercepting 
Surface - in/hr. 

A = ratio of evaporating surface to projectional 
area - sq.mi. 

ISC = storage capacity of the intercepting surface 


areas - in. 


SS = SURFACE STORAGE - in. 


DS = Depression Storage - in. 

DSC = Depression Storage Capacity - in. 

5 Koi refers to the infiltration Ge ae etic 
of the soils 

P = excess precipitation (P - IS) - in. 

B = ratio of depressional storage area to total 


drainage area. 


D = Detention (overland flow) Storage - in. 


L = length of the overland flow reach - ft. 
ics = infiltration characteristic of the soils 
q = rate of overland flow per unit area - in/hr. 
S = Slope of the overland flow surface - ft/ft. 
n = a Manning's roughness factor for overland flow 
1 = the inflow (precipitation) rate - in/hr. 
OLF = overland flow - in. 
I = Infiltration Rate - in/hr. 
SM = soil moisture -in. 
I, = initial infiltration rate - in/hr. 
a = infiltration capacity - in/hr. 
P = amount of precipitation excess - in. 
1. = inflow (precipitation) rate - in/hr. 
- Flow Equation - 
i = inflow rate - in/hr. 
K = hydraulic conductivity - in/hr. 
D = K/ d SM is the soil diffusivity - in/hr., where 
dv 
d (SM)= moisture capacity i.e. the differential changein soil 


dw moisture with a change in capillary potential 
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partial differential change in soil moisture 
with soil depth relative to the surface 
(when moisture supply is not a limiting 
condition) 


Saturated hydraulic conductivity Keer 


AERATION ZONE STORAGE (UNSATURATED SOIL PROFILE) 


soil moisture - in. per unit depth 
Capillary tension - in. 

gravity potential - in. 

total porosity of given soil profile 


unsaturated hydraulic conductivity. 


- Flow Equation - 


partial differential change in soil 
moisture with time, t 


partial differential change in soil 

moisture tension 

f (W); relates soil moisture as a function 

of the capillary tension w and 

f (W,SM); relates the unsaturated conductivity 
as a function of the capillary tension and soil 
moisture; these are secondary relationships 
required to be developed as aid in the numerical 
solution of the partial differential flow 
equation. 


INTERFLOW (INTERMEDIATE) ZONE STORAGE 


10 :Q 


the outflow rate from the IZS - in/hr., cfs 
relates to the volume of interflow in 
storage - cu.ft, 

total hydraulic potential. 


- Flow Equation - 
the Laplace's equation for unconfined steady- 


state flow in an isotropic porous medium 
saturated hydraulic conductivity. 


SATURATED ZONE STORAGE 


outflow rate from SZS - in/hr., cfs 

volume in storage in baseflow - in. 

total hydraulic potential (P/y + 2), where 
P is the pressure potential, y the density 
of water and Z, gravity potential. 
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SRO = SURFACE RUNOFF ROUTED INTO STREAMFLOW 


ior Surface discharge (output) from the basin 
at time Loy, in cfs 

X = Successive input volume up to time, t 

ep = a time-area distribution parameter (distribution 
graph) for the successive input volumes - e. g. 
I.U.H. technique; t = summation over time t, 
y 
a=o 


for increments, a. 


INFL = INTERFLOW ROUTED OUT OF IZS INTO STREAMFLOW 


St Oiag interflow rate - cfs at time, t 
K = an interflow recession constant 
i initial interflow rate (N. B. - routing 
relationship assumed to be similar to that 
of BFL). 


BFL = BASEFLOW ROUTED OUT OF SZS INTO STREAMFLOW 


dy a baseflow rate - cfs at time, t 
K = a baseflow recession constant - 

(N. B. - the routing function may require 
adjustments to account for variable 
component of ground-water recession) 

a. = initial baseflow rate. 


TRO = TOTAL OUTFLOW AT BASIN OUTLET 
- Route via linear channel storage and flow translation 
input 


output 
storage coefficient 


aL 
Q 
K 
dQ ‘ i 
ae time rate of change 1nQ. 


ET = ACTUAL EVAPOTRANSPIRATION 


SSM = supply of moisture (water) available for 
evapotranspiration - amount available in 
TSi'"'ss* “TZS. and: AZS = rin: 

= radiation = se 6 

= temperature - F, or C 

= wind- m.p.h. 

= relative humidity, H 

= a crop (vegetation) -use factor. 
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PRECIPITATION 


-ps Rain and/or 
Snow, P 
INTERCEPTION Interception Storage 
STORAGE- — --~-~ ~ Function 
- Is - 
Is = £(P, bi 1D Ne dete, 
Depression Storage 
Function 
£(Piy eC," 5) DSc) 
SURFACE Detention Storage 
STORAGE Function 
ao ie Dis of, ee, a) 
gqi= £ (L575) pe) 
Infiltration Rate Function 
Tee fiSM oP! feds, tele) 
e (0) Jel 
- Flow Equation - 
V =K+D _0(SM) 
o pam 
dy 
Soil Moisture Function 
(unsaturated flow) 
= K ia fds 
AERATION Pee War Seer o) 
ZONE { - Flow Equation - 
STORAGE cae oan @w + 1) 
- AZS - dy 
= £(W,SM) 
Interflow Function 
(saturated flow) 
INTERFLOW q= f(Q., H) 
ZONE : - Flow Equation - 
STORAGE 
=- IZS - V(K VH) = 0 
K =const., SM = const. 
sat. 
9 (SM) /dt = O 
Baseflow Function 
(saturated flow) 
SATURATED q = £(Q., H) 
ZONE ° 
STORAGE - Flow Equation - 
Sesys— V(K VH) = 0 
K = const., SM = const. 
sat. 
9 (SM) / % = 0 


Figure 4. HYDROLOGIC MODEL: 
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A Schematization for an Integrated (Comprehensive) 
Hydrologic System Model 
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HYDROLOGIC MODEL: A Specific Model Approach 
A Simple Regression Model 


The following gives an outline for the definition of a 
parametric, linear regression model, according to the conceptual 
model scheme of Figure 5a. 


_ The basic relationship for a simple linear regression, 
rainfall-runoff model may be expressed as: 


R= a+ bP secee te) 
; 2 i ran ee 
with b=. (Re i RD: oA Pe eee 4 2 (Pe =r) sec (LD) 
i=1 = Anette ge 
and a = R - bP 
where R = runoff 
P = precipitation 
a = least-square regression constant (usually 
negative) 
b = least-square regression coefficient (usually 
< tal 2:0) 
bh TS hoot I a a eae | 


Operation in the Whole System: One-Element Model EO - 
Figure 5a 


i) Relationship in terms of runoff, R: 


feo a” C 
re) 
eae is Saal 
inf 2 E> iC 
re) 
R= b a 
= oe 3 
R besBnch. a (39) 
with 
a =eb_C_ = constant. 
fo) ome) 


where, C_ is a numerical parameter in EO, related to the surface 


and subsurface storage or retention characteristic of the 


drainage basin, 
bo is regression coefficient, an operator needed to convert P 


(input to the system) to R (output from the system). and may be 


a a 


a 


V ~—2 


related to the runoff characteristics of the basin; a is a regression 


constant, or a parameter related to the moisture deficit and 
retention of the soils (a minimum moisture requirement of the 
basin to be satisfied before the initiation of additional runoff 
by a given input P). 


ii) Relationship in terms of basin loss, L: 


tae Pe eouwae 


= i] + i] 
L=b ae Peas ey) 
Weeden = (lL = b_) and 
O fe) 
' = = 
tats bo a constant. 
bts is a regression coefficient or a parameter needed to 


convert portion of the P (retained or stored in the system) to 
loss L (output from the system) and may be related to the 
transpiration and evaporation characteristics of the vegetation 
and soil, respectively, in the basin, 

aloe is a regression constant, or a parameter related to 
a minimum moisture loss from the system when the moisture supply 
isenot Limiting. 


Operation in the Whole System: Three-Element Model, 


FltjnE2 and E3 - Figure 5a 


i) In Element El: 


< 

If P - Cy 
= eee Gs 

Ry O (6a) 

= Tipe at 8 oh 5: 

Ly P (6b) 
cS ae ae 
1 

Rj = pp = Ci oawe Ue 

ite = (¢ e e (7b) 


paar’ ch at Beaige >. ee, tee) eee jaiblog F1tevnoS 
i 33 ig'¢? 6¢ Vat pos (eeteys off ek gouge) ee 
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- i. hace fa nt e: 
(3 = a0 
i) 7 (4) | 
j 3 7 


Ni ees 2S. 


where, Cy is a numerical parameter in Element, El, related to 


the initial surface storage or retention characteristic of the 
system, 


L. is the potential loss from Element, El, with a maximum 


1 
capacity C) and 
; Ry is the output from El which is a virtual input to Element, 
E e 


Element, El, may represent the surface storage characteristics of 
the basin (interception by vegetation is neglected). 


nip | In Element E2: 


OS b,R, exe (ia) 


ta =e. Db 


2 seaes COB) 


1) R 


it 
substituting equation 7a into equation 8a and 8b: 


Reb tb > C,) 


1 

R=b, P+ a, is ACEC) 

with 

a, = “b,c, = constant. 

Ly = La b,) Coes C,) 

Lead obs) 2 + a'y 

Ly = Daa t a's eau eC SC) 

with 

bY asp Oh aa b,) 

ay = “b') Cy = constant. 

where, by and by are factors, (regression coefficients) 
needed for the conversion of input Ry into runoff, R,or loss, 
Lo, respectively. 

ay and a's are regression constants or parameters which may 


be related to the moisture deficit and moisture retention that 
has to be satisfied, before the initiation of runoff or loss 


from the system. 
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Element, E2, may represent subsurface storage or retention 
characteristics of the basin. 


iii) In Element E3: 


L=.l, +L, hates toy) 
where, L is the total potential loss from the drainage basin, 
1. e. that portion of the precipitation input that is not accounted 
for by the measured runoff. 


E3 is a 'dummy' element introduced to provide for the sum- 
Mation of the potential losses from El and E2. 


Application of the Simple Regression Model 


Basic criteria: 


a) - drainage basin; select a basin that exhibits a recurrent 
storage condition at the onset of each major seasonal 
change, e. g. onset of major spring runoff when the 
basin storage is at a maximum; or select a basin with 
runoff sequences that may be regarded as having 
negligible contribution from adjacent basins. 


b) - time interval; select a time interval for the analysis 
that indicates a negligible lag or autocorrelation 
effect in the respective time series of precipitation 
or runoff data, e.g. monthly, seasonal or annual series. 


Data Requirement 


- concurrent series of precipitation and runoff(monthly or annual 
series). 

If monthly series are to be used, it is suggested that a 

normalization procedure be performed on the data prior to the 

regression analysis, e.g. use the logarithms of the original 

data assuming that the variates are log-normally distributed. 


Parameter Evaluation 


Gross estimates of Cy or Cy may be assigned from knowledge 


of the major soil types in the basin and their areal distribution 
and moisture characteristics, the surficial geological charac- 
teristics,the recharge and discharge areas and the major aquifer 
storage characteristics. 


The constants and coefficients in the regression model 
equation (la) may be determined by operating with the concurrent 
series of P and R in a standard least-square regression analysis 
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subject to satisfying a threshold precipitation value: 
R.=.a + bP, for ES ee 
R=o0, for P < Pp 
fe) 


where, Py is the threshold value (e.g. infiltration index) 


required to be satisfied before the initiation of Punocr, 


Expanded (Multiple) Regression Model 


Consider the multiple regression: 


where, y is the dependent variable (runoff) 

are the independent variables (precipitation, evapotrans- 
berotion, etc.) > K = 1,°2.2.en independent variables and 
a is the regression constant. 


If the basin water balance or conceptual storage equation is used 
as a physical base, 


yo) 
I 


ECP ily) 


with 1b 


fy (E, G) and 


Ss = f. (Si, S}) 


where, R is the runoff and is a function of the precipitation, P, 
the basin loss, L, and the basin storage, S; in turn, L is a 
function of the evapotranspiration, E, and the loss to ground- 
water flow systems, G; and S is a function of the soil moisture 
storage, Ss and the ground-water storage, Sp: 


Operation in the Whole System: Expanded Three-Flement 
miuecioie=rlement) Model’ — Figure 5b 


i) Operation in Elements El, E2 and E3 


It should be noted that the operations in El, E2 and E3 
are aS per previous development for the Three-Element Model, 
with the exception that the output Ro from E2 represent a virtual 


surface runoff. This is further augmented in E5 for baseflow 
contribution to give output R from the system: 
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Ro = biR, weed (20) 


where, the coefficient b, usually differs in value from the 


coefficient in equation (8a) of the Three-Element model. Also, 


the operations in E4a and E4b represent Sub-operations extruded 
from El and E2. 


ii) Operation in Flement F4 


By definition, S =L, that is, input to storage is equivalent 
to the potential losses, accounting for evapotranspiration, change 
in storages and loss to regional flow system. 


Let C., be a numerical parameter in Element E4 related to the 
total water-holding Capacity of the, soil sprotile, —.en) co bal 
porosity, and assuming that the surface supply is an instantaneous 
input to the soil profile. 


if Ss 5C5 
nn corr. (bd a) 
S, = 0 bee id 
Lies >Co 
Ss. a Ss. (capacity) Riexeeetlel G) 
Sp =S- Ss. (capacity) 


where, a is the increase to available soil moisture in the soil 


profile (unsaturated zone) and 
Sh is the percolation (drainage) to the water table 


(saturated zone). 
iii) Operation in Sub-Element, E4a 


By definition Ss. (capacity) = Car where, C3 is a numerical 


parameter in E4a related to the moisture capacity of the 
unsaturated zone, after drainage by gravity (i.e. field capacity). 


os > oF 

Past b.S wares Za) 
2a 

AS. = (l - bo) S5 Sie le) 


where, E is evapotranspiration (actual) — 
bo is a parameter or regression coefficient needed to 
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convert a portion of the available soil moisture to actual E, 
(1-b.) 1s related to the soil moisture-retention character- 
istics and 
AS, is the incremental change in soil moisture storage. 


iv) Operation in Sub-Element, E4b 


Let C, be a numerical parameter in E4b related to a base 


storage condition of the local ground-water flow system, i.e. 
the aquifer storage condition necessary to sustain a minimum 
baseflow (assuming that fixed proportions of the recorded 
precipitation are transformed simultaneously to baseflow and 
changesin storage). 


Le Sere Cy 
CC Se poVoa) 
G, = b36 eee oO) 
Ceszr° (assumed) Svat SC) 
AS) = (les b) G bie te a de oG) 
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G= opt tC Sere ta) 
Gh = b,G 
G = b,8, + a, Seta heb) 
with a3 = bCy = constant 
Genz © (assumed) 
and AS) = (a= b,)G 
bf oer Ce b) Sy + a', vated 4a) 
with a', Beth ir b3) CF constant. 
1s S774 
G = Sh at C, 
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with a". = bi 4c, = constant 
and AS) = (1 - db, = b',) G 

AS) = (1 - b, - b',) Ss, t+ at", are bated, 
with a'". = (1 - es - b' 3) C, = constant, 


where, G is the total potential outflow from the local ground- 
water flow systen, 
G, is the local ground-water contribution as baseflow, 


G. is the loss from the local ground-water flow system to the 
regional flow system through deep percolation, 


b, is the regression coefficient needed to convert a portion of 
the ground-water storage to a baseflow contribution and 
Cle b3) may be related to the aquifer storage characteristics. 


is a regression constant and may be related to minimum baseflow 
CONECELD ution, 


b', is a regression coefficient needed to account for that portion 

of recharge lost, through deep percolation, to the regional 

flow system and 

(l - b = b',) may also be related to the aquifer storage 

characteristics, 

a". may be related to a minimum deep percolation loss where 
recharge to local’ ground-water flow as not limiting, 

a'., and a'". may also be related to the overall aquifer storage 

characteristics, and 

AS, is the incremental change in ground-water storage. 


v) Operation in Element, E5 


RED RGB eG sForemt ELS.) 
where ,R igs the total runoff from the system. E5 is a 'dumny' 
element introduced to provide for the summation of the conceptual 
surface runoff and baseflow to produce the streamflow. 


Application of the Multiple Regression Model 


Consider the multiple regression model involving P, E and S 
as predictors of runoff, R: 


R= ay + bP + boE + b,8 Pa Lee 


where, aye Dy P, bo, Eg, b and S are as previously defined. 
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Definition of the Model Parameters 
IE ote tL Parameters 


Equation (16) can be written according to a regression model 
(Amorocho and Orlob, 1961) as equations (147 re 


m 
R=b 
By asP. Tice eee C157,) 


+ 
1 Py + by Et 5? 


2 5 
where, Ris the observed runor ie 
Py 1s the recorded precipitation, 


Eo is the computed estimate of evapotranspiration, 


Ee represents the storage term, S, in which P. is 
concurrent and antecedent precipitation 
VolLues 2e = ol 2)..44. period. 
ays cet bie bor are regression constants and coefficients. 
The input parameters P and EQ are determined with a certain amount 
of error in their estimates viz., 
a) Inflow - Index: 


— 1 
PL b 1? + en wee 1OLS) 


where, P, is an estimate of the areal average precipitation, 


P is the station(s) recorded amount, 
by a station weight (coefficient) for areal distribution of P 


and 
e> is the standard error of the average. 
b) Evapotranspiration - Estimate: 
= J ate 
EQ b 9E + en sp he MEI 


where, E_ is the areal average estimate of evapotranspiration, 
E is the station(s) computed estimate, 
by a station weight (coefficient) for areal distribution of 
E and - the standard error of estimate. 


c) Storage - Index: 


If it is assumed that fixed proportions of the total recorded 
precipitation are transformed into ground-water outflow, sim- 
ultaneously producing changes in storage, then the storage index 
term of equation (17) may be interpreted as: 


= ee 20 
+ G. ae ay (20a) 


where, AS_, AS, are change in soil moisture and ground-water 
. storage and GC. is deep percolation loss. 


a b 


If G. is assumed to be negligible, with an appropriate selection 
of the time interval, the storage index, S, may be estimated 
with presumably small error, from: 
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Sey 7 & eo (OUTS) 


where, In is the net inflow at the soil surface and G, the 
baseflow, 
which may be determined from total runoff hydrographs by hydrograph 
separation techniques. Incidentally, the net inflow may 
be examined on the basis of: 
m 
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Model Simplification 


_The following statistical criteria are given for the sim- 
plification of the model relative to the selection of appropriate 
time periods for the analyses. At certain periods, some of the 
major parameters may tend to a constant due to recurring climatic 
conditions. For example, to examine the statistical degree of 
accuracy with which R can be estimated on the basis of a 
multiple regression model, 

n 
‘Lage zr 1 xy +a era) Cee eh) 


Where, X, are random independent variables), 


usually the best choice of an end-point is when the variance of 
the observed values of a given variable is at a minimum, or when 
its values are approximately equal to the average value for the 
period of observation. 


To examine this, consider the variance of the regressed 
estimate, y: 


War AV). = Var AX.) 2 Fin COM Jl Mea te Set) 
k ; 3,k ues 
where, X. and Xy are specified independent variables. Minimization 


of the error in y may be examined by the following relationships: 


0 Var ( ) ~ (23) 
Var (x) racer on 

09 Var (y) =O Stee Cee) 
0 Cov (XX) 


thateisgntaking the partial’ derivative of the variance of y with 


respect to the variance of X and the covariance of Rs and Xi, 


and equating to zero, respectively. 
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As the sample estimate of Var (X)) is given by: 


Wat (xa) e ae sf cep cy 
k eee (25) 
Nie 4b yk us : 
eye ca ge 
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Hence X, a x for minimum variance satisfying the partial 


derivative given by equation (23). Similarly, the sample 
estimate of the covariance is given by: 


COV aX: X = 1 / 
5 k? y (XX 
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Hence A > xX, for minimum variance satisfying the partial 


derivative given by equation (24). 


These criteria may be applied to model equation (16). 
Consider a time interval such that the values of two of the 
independent variables approach their respective mean values, 
for each set of observations, 


Viz. E> 'E 
c Cc 


AS + AS 


or choose a period when the mean value of one variable becomes 
independent of the magnitude of the other variable, i.e. the 
conditional expectations of one variable's values, given the 
other variable, tends to a constant: 


E (E | AS) >. constant + E. 


(ong 


E (AS|E,) + constant > AS 
A similar consideration may be given to other pairs of 
independent variables. 


By operating with the above criteria, the relationship 
(equation (16)) may be simplified: 


(a) for an observational period with a constant evapotrans- 
piration rate, 
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where, CE is a constant evapotranspiration rate for At, the time 
interval (ty - ty): 


(b) for an observational period with approximately equal changes 
in storage, - considering each storage component separately - 


2.) for baseflow: 


AS, >AS, > Oo; Var (AS,) > fe) 


and G.> G SyC ht Var (G) > fe) She. 1030) 
where, CL is a constant baseflow rate for At, the time interval 


ii) for a period with constant AS. (period when near maximum 
saturation or near complete depletion of moisture in the 
aeration zone can be assumed): 


AS, 7 AS, > CO. iV ar (AS_) > re) 


Now considering together conditions (a) and (b), the 
equation for surface runoff R, may be expressed: 


Roe b Py = (C, + C,) At Hay ces Aes) 
Le BE. and G, cannot be estimated independently, then a 
general relationship between runoff and precipitation 
would result with regression coefficients that reflect 
the effect of evapotranspiration and baseflow. 


In general, on a water-year basis with a carefully 
selected end-point, e. g. latter part of the drought 
season, when 


AS, > AS), 


Wie oo. Oo; Var (AS_)+0 
a a a 


> oO; Var (AS) > fe) 


and 


GC. = o (assumed) 


the relationship may be expressed as: 


és a . Sag CST 
R b, Py Ee (32) 
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1 6 Ea is not estimated independently; the regression is: 


R= b,P, + a, pive (33) 


in which case a regression would be the result with the effect 


of Eq reflected in ay and by 


Data Requirement 


- historical record of the various variables. 
For an annual P-R model analysis (arranged on a yearly basis, 
beginning with each month), a concurrent series is required: 


Lee - 
4 Sb 
2-R 5 = 8 
a 
3 - Ey 6 - P (snowpack storage) 


For monthly P-R model analysis, the following concurrent series 
is required: 


ee ai Si, (at the beginning of the month) 
2 = R 6 - 5 (at the beginning of the month) 
3 = EQ 7 - Sy (at the end of the month) 
4 - P_ (snow) 8 - Ss (at the end of the month) 
95 other meteorological parameters 
(temperature, humidity, radiation). 
Parameter Evaluation 
Estimates of Cy Co and C. may be assigned from knowledge of 


the major soil types in the basin and their areal distribution and 
hydraulic and hydrologic characteristics; estimates of Cy may be 


assigned from knowledge of the basin surficial geological char- 
acteristics, recharge and discharge areas and the major aquifer 
hydraulic and hydrologic characteristics. 


The constants and the coefficients of the regression models 
(ewazeequations.31, 32, 33).may be determined by least-square 
techniques through operations in the appropriate regression 
Sub-routines. These regression sub-routines are the result of 
the appropriate application of the statistical criteria for the 
simplification of the multiple regression models equation 16 or 
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Whole System: One-Element Model 
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Figure 5a: 


Whole System: Three-Element Model 
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HYDROLOGIC MODEL: A Definition Scheme for a 
Parametric, Linear Regression Model 
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Figure 5b: 
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Whole System: Expanded Three-Element 
(Multiple Element) Model 


Symbols 


P - precipitation 
R, - net input 


R Ro runoff 


Lie Lo, L - loss term 

iS) = STOrage 

S. - soil moisture storage increment 
Sy), - ground-water storage increment 
G, - baseflow contribution 

Gc. - loss to regional flow 

E - evapotranspiration loss 
El,...,E5 - element of the system 


HYDROLOGIC MODEL: 


A Definition Scheme for an Expanded, 
Parametric, Linear Regression Model 
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Figure 5c: 


HYDROLOGIC MODEL: Flow Chart for Parametric, 
Linear Regression Models 


Whole System: One-Element Model, EO 
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INPUT DATA: 


PARAM., EQN. - 
Coef. & Const. 
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Vary 


Whole System: Three-Element Model, E1,E2 & E3 
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Whole System: Expanded Three-Element (Multiple-Element) 
Model, El, E2, E3, E4, (E4a, E4b), E5 


START 


PARAM. , EQN. 
Coef. & Const. 
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LIST OF SYMBOLS IN MODEL (FLOW CHART) 
ONE-ELEMENT MODEL EO 


Py = recorded precipitation - in. 

R, = recorded runoff - in. 

c. = parameter for basin retention - in. 
aor a' = regression constants 

bor De = regression coefficients 

i=l, 2, 3,....N periods of observation 
Ry = computed runoff - in. 

L, = computed loss - in. 

Li regression residual - in. 

ISE = integral square error 


THREE-ELEMENT MODEL El, E2, & E3 


Ps = recorded precipitation - in. 

R; = recorded runoff - in. 

Cy = parameter for basin retention - in. 
apr aly = regression constants 

by, b'y = regression coefficients 

Li; = a loss term, from El — in. 


computed net input to E2 - in. 
= a loss term, from E2 - in. 


R 
L 
L; = sum of the losses - in. 
Rie computed runoff -in. 

i. 


Eee, 3....N periods of observation 


MULTIPLE-ELEMENT (Expanded 3-Element) MODEL 
Element E4, F4a, E4b & E5 (see definitions above plus) 


S. = storage input (= potential total losses) - in. 
Co = a parameter for basin retention - in. 
Cc, = a parameter for basin retention - in, 
S4im soil moisture storage ingrenene as in; 
Shi = ground-water storage increment - in. 


es evapotranspiration - in. 
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AS 43 = change in soil moisture storage - in. 
Cy = a parameter for basin storage - in. 

G; = a potential outflow increment - in. 
Gos = a baseflow contribution - in. 

Gad = loss to regional flow system - in. 
ASLi = change in baseflow storage - in. 

Roy = a surface runoff component - in. 

R; = computed runoff - in. 

ass a's, anys am 4 = regression constants 

b be, b', = regression coefficients. 
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DROLOGIC MODEL: A Brief on Stochastic 
Models 


Figure 6: 
Stochastic Hydrologic System 


M Model 


Tom wi ui? ih 
oh | 
| a 
an 
+ car F 


€ 


Page 
VI-1 


A Schematization for a Type of 
VI-3 


4 


yt 


a eal 


ae i POF ee wae 
i: 7 ya dire. 


ne hedip “ | 
tee Pate 


ay Sis 


. 
Ne Ai. 


| 
| 
) 


HYDROLOGIC MODEL: 
A Brief on Stochastic Models 


The behaviour of a hydrologic system and i 
with a sequential time menor. this ceili ee 
degree of uncertainty which can be described by probability laws 
and, as such, the System or process may be regarded as stochastic 
(an index family of random variables). That is, there is an element 
of chance in the values which the variables may take in the 
interrelated hydrological processes which define the average 
behaviour of the system. 


Parametric hydrology generally deals with discrete hydrologic 
events, such as storms or floods. A parametric approach is 
deterministic in the sense that derived functional relationships 
between the input-output processes attempt to define the behaviour 
of the system for all times; that is, a given set of causes always 
produces the same set of effects. Stochastic hydrology, on the 
other hand, concerns itself with the time-sequential properties of 
the discrete events. It attempts to formulate a mathematical 
model for an observable phenomena which changes in time, in a way 
that is not completely predictable, but is expressible in terms 
of probabilities. 


A conventional approach to stochastic hydrologic system 
modelling usually takes the form of a hybrid model (deterministic 
and stochastic). To achieve this hybrid model, a firm foundation 
based on the deterministic model concepts is required. The model 
may be used to transform the recorded series input of a variable 
with known or assumed probability distribution to values of an 
output with unknown probability distribution. Time~series analysis 
may then be applied to examine the generated output for its 
stochastic properties. Alternatively, the stochastic properties 
of the input variable may first be determined, then a deterministic 
relationship is used to transform the input to the required 
output. The generated output is presumed to have stochastic 
properties similar to that of the input variable. 


Another alternative method, such as in the case where a 
deterministic relationship is known to exist between two time 
series, but cannot be specified explicitly, is to establish 
through regression or correlation, empirical relationships between 
parameters derived from statistical analyses of the individual 
time series. These empirical relationships are then used as a 
deterministic base, for the transformation procedures. 


A more aggressive approach (Chow & Kareliotis) to stochastic 
hydrologic system modelling, Figure 6, is to formulate a stochastic 
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Views 2 


model for the hydrologic system in which runoff is considered as 
the integral of three component stochastic processes, viz., 
change in conceptual watershed storage, the total rainfall input 
and the total watershed losses, mainly evapotranspiration. In 
the analysis, the integrated stochastic processes are treated 

as a three-dimensional vector or multiple-time series with each 
time series consisting of a deterministic and a random component. 


The mathematical notation given for each stochastic process 
in Figure 6 is called a time series. For example the set of 
observations (P (t); tBT) is called a time series for precip- 
dtation; i.e. for each t in T the observation P (t) is an 
observation of the random variable. Prior to analyzing the time 
series, it is necessary to assume a model for the stochastic 
process, e. g. moving average, MA, harmonic, H, auto-regression, 
AR, Or Markov Chain, MC. Standard statistical techniques in 
time-series analyses (correlogram and spectral density analysis) 
aid in the identification of the model inherent in the data 
sequence, through goodness of fit tests on the assumed model by 
comparing generated output with historical data sequences. 
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